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We present a new method to study light scattering on non-absorbing spherical 
particles. This method is based on the Ramsauer approach, a model known in 
atomic an nuclear physics. Its main advantage is its intuitive understanding 
of the underlying physics phenomena. We show that although the approxi- 
mations are numerous, the Ramsauer analytical solutions describe fairly well 
the scattering phase function and the total cross section. Then this model is 
applied to the Henyey-Greenstein parameterisation of scattering phase func- 
tion to give a relation between its asymmetry parameter and the mean particle 
size. 
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The subject of light scattering by smallparticles is present in several scientific areas such 
as astronomy, meteorology and biology [l, 2|]. The first model by Lord Rayleigh in 1871 dealt 
with light scattering by particles whose dimensions are small compared to the wavelength. 
A significant improvement came with the Mie solution {3] which describes light scattering 
by spherical particles of any size. In an apparently different domain, the scattering of a low 
energy electrons by atoms studied by Ramsauer experimentally [4] showed surprising struc- 
tures, and it was several years before the explanation found by Bohr describing the electron 
as a plane wave (5 T] ! The Ramsauer effect was certainly the first phenomenon showing 
the wave properties of matter. This framework has since then been used to describe very 
different types of collisions such as atom-atom [8j or even neutron- nucleus 9Nl3| . Our idea 
is that since light behaves like a wave, it should be also possible to apply Bohr/Ramsauer's 
ideas to light scattering by spherical particles. 

Sect. [1] is a brief introduction of some quantities concerning light scattering on (non- 
absorbing) spherical particles. Sect. |2] describes the Ramsauer effect and provides details on 
the way to obtain the Ramsauer scattering phase function and its simplified solution. We 
compare the Ramsauer predictions with the Mie calculations, emphasizing certain charac- 
teristics of light scattering. Then, in Sect. [31 we integrate the differential cross section to 
obtain the total cross section and the associated extinction efficiency parameter using the 
Ramsauer approach. Finally, Sect. H] presents an application of the Ramsauer approach for 
the Henyey-Greenstein parameterisation of the scattering phase function; a basic relation 
between its asymmetry parameter and the mean particle size is obtained. 

1. The Mie predictions for scattering of light over a sphere of non-absorbing 
dielectric 

When light interacts with particles, two different kinds of processes can occur. The energy 
received can be reradiated by the particle at the same wavelength. Usually, the reradiation 
takes place with different intensities in different directions. This process is called scattering. 
Alternatively, the radiant energy can be transformed into other wavelengths or other forms 
of energy, such as heat: absorption takes place. In the following we will consider absorption 
equal to zero. The intensity of light scattered at an angle 9 and a distance r by a single 
spherical particle with a radius R, from a beam of unpolarised light of wavelength A and 
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intensity Iq is given by 



I{6) = ^xa sca (R,\)xP(d), (1) 



r 



where Jo is the initial intensity of the beam, a sca (R, A) the scattering cross section and P{9) 
the scattering phase function. The integral of P(9) over all solid angles has to be equal to 
unity. The scattering cross section a scai (R, A) has the dimension of area but is not in general 
equal to the particle cross-sectional area. Indeed it is customary to define a scattering 
efficiency 

_ o~ sca ,{R, A) 

where ttR 2 is the geometric cross section in the case of a sphere of radius R. The scattering 
efficiency represents the ratio of the energy scattered by the particle to the total energy in the 
incident beam intercepted by the geometric cross section of the particle. Since absorption 
is assumed null here, the extinction efficiency Q ext is equal to Q sca and the corresponding 
total cross section <j tot is equal to <r sca . 

As already mentioned, the ratio of the particle size to the wavelength of the radiation of 
interest is of crucial importance for the particle's optical properties as well as for the choice 
of a suitable method for calculating those properties. We therefore define the dimensionless 
size parameter as 

2ttR , 

X = —j— = fcout R, (3) 

where R is the radius of a spherical particle, A the light wavelength and k ont the wave number 
outside the sphere. In the case of non-spherical particles, R might represent the radius of a 
sphere having the same volume. The theory describing the scattering of an electromagnetic 
plane wave by a homogeneous sphere of arbitrary size was originally presented by Gustav 
Mie [3]. Particles with a size comparable to the wavelength of visible light are relatively 
common in nature. In this case, the Rayleigh theory is not applicable anymore because the 
field is not uniform over the entire particle volume. The laws describing the total scattered 
intensity as a function of the incident wavelength and the characteristics of the particle are 
much more complex than those for Rayleigh scattering. 

A. Extinction efficiency: the total cross section 

The 1/A 4 dependence of the total scattered intensity found by Rayleigh is no longer valid 
generally for the Mie solution. The extinction efficiency factor Q ext = a tot /(iiR 2 ) depends 
on the radius R and on the relative refractive index of the particle n. The relative refractive 
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index is denned as the ratio of the refractive index of the particle over the refractive index of 
the medium. In the following, where the medium is air, the relative refractive index corre- 
sponds to the refractive index of the particle n. The Mie theory uses Maxwell's equations to 
obtain a wave propagation equation for the electromagnetic radiation in a three dimensional 
space, with appropriate boundary conditions at the surface of the spherical particle. The 
extinction efficiency factor obtained is 

2 00 

Q cxt = ^^(2£+l)Re(a, + 6 £ ), (4) 
i=\ 

where the Mie scattering coefficients and bp are functions of the size parameter x = 
2txR/\ = k out R and of the index of refraction n jl, 2]. To compute Q ext numerically using 
Eq. (j4j), it is necessary to truncate the series, keeping enough terms to obtain a sufficiently 
accurate approximation. The criterion developed by Bohren in (2( was obtained through 
extensive computation. Analysis of the convergence behaviour of Eq. (T4J) revealed that the 
number of required terms N was slightly influenced by refractive index and had to be the 



closest integer to x + 4a; 1//3 + 2 14J. For instance, for a raindrop of 50 fim radius and a visible 
wavelength of 0.6 /im, the number of required terms is N = 558. Nowadays, computers 
have reached a point where computing time is no longer a major problem. The closed form 
Ramsauer approach, which we will develop later, would have been useful some years ago 
to provide computational ease. In the present context it furnishes useful insights into the 
different physical processes involved. 

Fig. [1] shows, for three values of the index of refraction, Q ext as a function of the size 
parameter x. Each curve is characterised by a succession of maxima and minima with 
superimposed ripples. The amplitude of the large oscillations and of the ripples grows with 
n. For particles much larger than the wavelength of the incoming light, i.e. with a size 
parameter x ^> 1, it might be expected from geometrical optics that the total cross section 
would be irR 2 , and thus Q e xt would approach 1. However, Fig. [1] shows an interesting 
and somewhat puzzling trend, the scattering efficiency approaching 2. It is the so-called 
extinction paradox, and is discussed in detail by Van de Hust The phenomenon is 
apparent only for observations made far from an object, so that light that is scattered 
even at a small angle can be considered removed from the beam. The total cross section 
results from the geometrical contribution (classical limit) of nR 2 and an equal contribution 
from the diffraction of the incident plane wave at the surface of the sphere. The diffraction 




Size parameter 

Fig. 1. (Color online) The extinction efficiency factor Q ext versus the size parameter x, for 
non-absorbing spherical particles with relative refractive indices n = 1.05, 1.33, and 2.00. x 
is given by the relation x = 2ttR/X = k out R. The vertical scale applies only to the lowest 
curve, the others being successively shifted upward by 2. 
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Fig. 2. (Color online) Angular scattering for water droplets illuminated by unpolarised 
light, for three different radii. The incident wavelength is fixed at 0.6 /xm. The vertical scale 
applies only to the lowest curve, the others being successively multiplied by 100. 
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Fig. 3. (Color online) Angular scattering for a distribution of water droplets illuminated 
by unpolarised light, for three different standard deviations. The mean radius R is fixed at 
15 //m and the incident wavelength at 0.6 /zm. The vertical scale applies only to the lowest 
curve, the others being successively multiplied by 100. 

contribution is strongly peaked forward. For nearby macroscopic objects, this diffracted light 
is not distinguishable from unscattered light at 9 = 0° and the paradox is not observed. 
B. Angular scattering: the differential cross section 

Mie scattering by single particles irradiated by laser sources is sufficiently strong to be 
detected with high signal-to-noise ratios for particles larger than about 0.1 /zm. The noise 
results from Rayleigh scattering by residual molecules and from the electronics used. The 
recorded signal depends on the scattering angle, as well as on the particle size and refractive 



index 



15]. The angular dependence can be calculated from Mie theory. In this study, we 



use a code developed by Philip Laven |l6|. For values of the size parameter x approaching 
unity, an asymmetry favouring forward scattering appears. For j;> 1, forward scattering 
increases even more strongly, showing very rapid changes for small increases in the scattering 
angle 9. Some of this behaviour is shown in Fig. [2] for water droplets of different sizes when 
illuminated by unpolarised light of A = 0.6 /im. Very small droplets show typical Rayleigh 
behaviour with a symmetric curve and a weak minimum at 90°. When the size parameter 
x increases, additional minima and maxima appear and a strong asymmetry develops with 
the forward scattering several orders of magnitude stronger than the backscattering. 
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Most natural and artificial collections of spherical particles do not exhibit the large fluc- 
tuations as in Fig. |2j Even a narrow dispersion in the distribution of the particle size washes 
out these features that strongly depend on size. The size distribution can be described using 
a log-normal function defined as follows 

n(R\R, a) = = -=- exp — — — 2 , (5) 

dit V27r log 10 a R V 2 log 10 a J 

where R is the radius of a spherical particle, a is the geometric standard deviation, R is the 
geometric mean radius, and N is the total number of particles. Phase functions for three 
different geometric standard deviations are plotted in Fig. [31 with a mean radius R fixed 
at 15 /im and a wavelength equal to 0.6 fim. As expected, a smaller geometric standard 
deviation implies a wavier scattering phase function. The ripples are stronger in the forward 
and backward scattering peaks. Consequently, in nature, where particle radii are always 
dispersed, light scattering will be described using smoothed phase functions (see Sect. HJ). 

2. Application of the Ramsauer effect for light scattering 

The Ramsauer effect was discovered in 1921 J4] while studying electron scattering on argon 
atoms. The total cross section versus the electron energy showed a surprising dip around 
1 eV. In Fig. H] we see recent measurements of electron scattering over krypton and neutron 
scattering on lead nuclei. 
A. The Ramsauer scattering phase function 

The idea in modelling the phenomenon is to consider the incident particle as a plane wave 
with one part going through the target and another which does not (Fig. EJ. The two parts 
recombine behind the target and then interfere with each other, producing the oscillating 
behaviour. Depending upon the impact parameter b (Fig. EJ), light rays going through the 
drop accumulate a phase shift. 

Using the Huygens-Fresnel principle, the propagation of light can be described with the 
help of virtual secondary sources. Each point of a chosen wavefront becomes a secondary 
source and the light amplitude at any chosen location is the integral of all these sources 
emitting spherical waves towards that observation point. The secondary sources are driven 
by the incident wave. Their amplitude is 

i (1 + cosfl) 

— ^ AncidontCLO [O) 
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Fig. 4. (a) Electron- krypton normalised total cross section versus energy [7|. (b) Normalised 



neutron-lead total scattering cross section versus energy 
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where A is the wavelength, Ancident ^ s the amplitude of the incoming light, 9 is the angle 
between the direction of the incident light and the direction from the virtual source to the 
observation point, and dS is the elementary area of the secondary source stands. 

We define the plane P, shown in Fig. [6j as the location of our secondary sources. The 
shadow of the sphere on the plane P is the disk CD. If the plane Q is chosen as the origin 
of the phases, the incident light is a plane wave which, on the plane P, has the amplitude 
e l2koutR . The amplitude distributions of the secondary sources can be split into two parts: an 
undisturbed plane wave and a perturbation in the CD region only. The undisturbed initial 
plane wave, when integrated, is described by a Dirac delta function in the forward direction. 
Thus the differential cross section comes from the sources on CD. The observation points 
are considered to be very far away so, to get the amplitude at an angle 9, we just sum all 
directions parallel to 9. 

Let p and (f) be the polar coordinates in the plane P. The amplitudes from the virtual 
sources in the disk CD are then 

1 + cos 9 



Aiisk(p) 



i 

A 



J2k out R , e in2k out ^R 2 -p 2 e «A: out (2R—2y/ R 2 -p 2 ) 



■ ^out i2k out R 

2tt 



i2k outy /R 2 -p 2 (n-1) 



1 + cos 9 



(7) 



The resulting scattering amplitude /sphere (#) is obtained after moving the origin of the 



8 




(a) 



Fig. 5. (Color online) Qualitative picture of the Ramsauer phenomenon. The wavelength of 
the light is supposed to be reduced in the dielectric. This picture shows the case where the 
contraction of the wavelength between (a) outside and (b) inside the medium (dashed line) 
is such that they come out in phase. Thus the sphere looks invisible resulting in an almost 
zero cross section. 

phases to the plane P. This is simply done through a multiplication by e~ l2k ° ntR 

-2tt r R 



/sphered) = e^*"** / d0 / A disk ( P ) p ^dp 



hut 1 + cos 9 



JO 

2tt 



d<p I < 
i) Jo 



R 



jfcoutpcos <fi sin t 



1-e 



i2k out ^R 2 -p 2 (n-1) 



(8) 



pdp, 



2tt 2 

where the term pcos0sin# represents the path length difference between a source in the 
disk and a source at the disk centre for an observer placed at an infinite distance from the 
disk. 



The integral over the angle is a Bessel function 



Thus Eq. (jSj) becomes 

R 

Y _ e i2k outy /R 2 -p 2 (n-1) 





r fn\ -1 1 + COS^ , 
/sphere (0) = «&out ~ / P M&outP Sin ( 



(9) 



dp. 



This integral has no analytical solution. The Ramsauer differential cross section is 
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Fig. 6. Definition of the variables used in the text to calculate the phase function for light 
scattering by a spherical particle of radius R. k in and k out are the wavenumbers for the light 
inside and outside, respectively, and b is the impact parameter. 



and we deduce directly the Ramsauer phase function 
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a 



da 
d9 



I /sphere (0)1" 



2n 7r 



(12) 



sphere fd$fd9sm6\f sphc 





This solution has to be computed numerically. Thus, this illustrates that predicting the 
behaviour of light scattering is not easy. 

B. The simplified Ramsauer scattering phase function 

We can get a formula in closed form if the sphere is approximated by a cylinder of radius R 
and of height 2R along z. Eq. ( TTUj) simplifies into 

l + cos# fR 



/cyl(0) = ^out" 



1 + COS0 
'"'out L 1 



/ pJo(Wsinfl) [1 - e «*out*(«-i)] dp 
Jo 

f R 

] / pJ (^outpsin6')dp 
Jo 



J2{n-l)k ont R-\ 



(13) 



.1+COS6 1 , 2(»-i)fc out a . 

i e 2 (— 2z)sm 



2(n - 1) k out R 



R 
sin 9 



Ji(k out Rsmt 



since the integral can be simplified using the fact that 

R 



J uJ {u)du = RJ 1 (R). 



(14) 
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The differential cross section is 



dcr 
d9 



I/. 



cyl 



R 2 



cyl 



1 + cos^ m 2J,{k mt R sin I 

^out -fl, ^ Sill [^J4 1 JAtoutJXJ 



(15) 



2 1 ""' J /.-,, ,: /)' -.in 9 

In fact, in the Ramsauer approach, the total cross section for a cylinder of radius R and of 
height L is given by 



2tt 



a 



cyl 
tot 



d$ / d#sin# 



da 
d9 



AttR 2 sin 2 



cyl 



fn-r 



AnR 2 sin 2 [(n - l)fc ou t#] 



where L = 2R in our case. Thus, the differential cross section can be expressed as 



da 
d9 



cyl 



1 

47T 



a 



cyl 
tot 



knntR' 



1 + cos#) 2J 1 (fc out i?sini 



(16) 



(17) 



2 k out Rsm9 

As we will see in Sect. 12 C\ it is possible to use a similar equation for the differential cross 
section with a sphere, i.e. 



da 
d9 



sphere 



1 

47T 



a 



sphere 
tot 



kontR' 



1 + cos9) 2Ji(fc out -Rsin( 



2 k out Rsin9 " <N) 

In this subsection, we have shown that in the simplified Ramsauer approach (SR), the 
total cross sections for a cylinder and for a sphere can be different but the scattering phase 
functions coincide and are given by the relation 



P (9) - 1 da - 1 



k utR 



1 + cos 9) 2J 1 (fc out i?sini 



1 2 



(19) 



2 k out R sin 9 

It seems interesting to note that this simplified phase function is independent of the relative 
refractive index. The validity of this independence will be verified in Sect. 12 CI It has been 
assumed in this last part a cylinder of radius R and of height L = 2R, i.e. a volume equal 
to 2irR 3 . But, in the case of a sphere, the volume of the scattering center is different and 
equal to (4/3) 7rR 3 . However, since the differential cross section is normalised to the total 
cross section to get the phase function, this difference is avoided. 

As can be seen in Eq. (Tl9"]) . the Ramsauer approach is much more intuitive than Mie 
model. We note two properties of the Ramsauer phase function: 

1. The amount of light scattered in the forward direction (9 = 0°) is proportional to the 
target area, i.e. square of the target radius: at constant wavelength, a larger sphere 
scatters more light in the forward direction, in agreement with the Mie prediction. 

1 dcr 



crdtf 



(9 = 0°) oc R 2 . 
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2. The full width at half maximum (FWHM) of the forward scattering peak is propor- 
tional to the inverse of k ou tR. 

2 



AO 



FWHM 



By detailed Mie calculations, these phase function properties deduced from the analytical 
Ramsauer solution can be verified. Fig. Ogives the phase function for pure forward scattering 
as a function of R and the full width at half maximum as a function of the size parameter. For 
three different refractive indices, the observed behaviour is confirmed by the Mie calculations. 
However, an analysis comparing the mean radii estimated by the Ramsauer approach with 
the true values used in Mie calculations has stressed an overestimation of around 30% of the 
values. This overestimation is independent of the incident wavelength and of the relative 
refractive index. After all, whereas it is difficult to obtain clear physical interpretations 
with Mie theory, the solutions being in the form of infinite series, the simplified Ramsauer 
solution makes easier prediction of light scattering behaviour possible. 
C. Validity of the Ramsauer solutions 

From the Ramsauer approach, two scattering phase functions were obtained: one from an 
exact calculation, Pr(9), and the second resulting from an approximation but given by a 
closed formula, -Psr(#)- The goal of this section is to evaluate the two solutions with respect 
to the Mie theory, taken here as the reference. As noted above, the simplified solution 
does not depend on the refractive index of the target. It is thus important to see how 
the Mie calculations and the Ramsauer solution evolve with respect to the refractive index. 
Fig. [8] shows the phase functions obtained from (a) the Ramsauer solution and (b) Mie 
theory for different refractive indices. For comparison, the one obtained from the simplified 
Ramsauer formula is also plotted with a thick black line. Although the phase function is 
almost independent of the refractive index for high n values, this is not the case when \n — 1| 
approaches 0. Moreover, the Ramsauer solution is in better agreement with the Mie curves 
when the refractive index is low. 

Consequently, when we want to study light scattering phenomena for spherical particles, 
as a first approximation, we can develop models without knowing the refractive index of the 
particles. This is due to the fact that Fraunhofer diffraction is mostly responsible for small- 
angle scattering and it is nearly independent of the relative refractive index. Thus, forward 
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Fig. 7. (Color online) Characteristic behaviours of physics quantities from the phase func- 
tion. The curves are obtained by Mie calculations (using [16]) for an incident wavelength of 
0.6 /nil. Three different refractive indices are shown: n = 1.05 (blue circles), n = 1.33 (red 
squares) and n = 1.60 (magenta triangles). The vertical scale applies only to the lowest 
curve, the others being successively multiplied by 100. 
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Fig. 8. (Color online) The simplified Ramsauer solution, independent of the refractive index, 
compared to the Ramsauer solution and Mie calculations, (a) SR vs Ramsauer. (b) SR vs 
Mie theory. Three different refractive indices are shown: n = 1.05, 1.33 and 1.60. The 
simplified solution, independent of the refractive index, is in black thick line. The incident 
wavelength is fixed at 0.6 fim and the sphere radius is equal to 10 fim. The vertical scale 
applies only to the lowest curve, the others being successively multiplied by 100. 



14 



scattering is sometimes favoured in the design of optical particle counters to eliminate the 
effect of refractive index on the measurement of particle size. 

The different expectations deduced from the Ramsauer approach could be also checked by 
a basic experimental setup. One of the instruments for angular light scattering measurements 



is the polar nephelometer [17|, [l8j. It consists of a collimated light source (usually a laser) 



and an arm that can be rotated around the sample of scattering particles. Two kinds of 
light collectors could be used: a single detector adjustable to different angles, or a set of 
multiple fixed detectors positioned around the sample and collecting simultaneously the 
light scattered at different angles. If the main goal is to get the best angular resolution for 
the scattering phase function, the first solution would be employed. The sample could be 
composed of polystyrene beads in suspension in water. Different mean sizes and different 
size dispersions are available in 19j. We note that a detector linearity is required over 



a large intensity range. Indeed, as shown in Fig. |2j a high intensity difference can be 
probed between the forward scattering peak and the rest of the angular range. Using the 
experimental setup briefly presented here, it seems possible to check the different theoretical 
expectations deduced by the Ramsauer approach. However, this work is beyond the scope 
of this study but should be done in future. 

3. The Ramsauer total cross section and extinction efficiency 

Since the results of the phase function predicted by the Ramsauer approach seem in agree- 
ment with the Mie calculations, this section will check if this is also the case for the total 
cross section. Two different methods will be presented to get the cross section by the Ram- 
sauer approach: the first one based on the optical theorem (and independent of the results 
given in Sect. [2]), and the second by integrating the differential cross section over (6, 0). 
A. Ramsauer total cross section via the Optical theorem 

In scattering theory, the wave function far away from the scattering region must have up to 
the normalisation factor the form of a plane wave and an out-going spherical wave, 

*(f) = e ik °" tF +f(0) . (20) 



The optical theorem asserts that 



4:71 

a tot = —lmf(9 = 0), (21) 
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where f(6 = 0) is the forward scattering amplitude. Under the approximation for the 
scattering of a scalar (spinless) wave on a spherical and symmetric potential, the scattering 
amplitude at a given polar angle 6 can be written as a sum over partial wave amplitudes, 
each of different angular momentum £, 

1 oo 

f W = VT~ 5> £ + he™' - 1] , 



(22) 



out 



£=0 



where rjt is the inelasticity factor (rji = 1 in our case of a non-absorbing sphere), Si is the 
phase shift (Si is real for pure elastic scattering), and P^(cos#) is the Legendre polynomial. 

From Fig. |6j St = (ki n — k ont )Rcosip. In the approximation of forward scattering, 
Pe(cos8) —7- 1. From Bohr momentum quantisation applied to the impact parameter 
b = Rsinip (see Fig. [6]), we have bp = £h and p = hk ont with £ = bk ont . By substituting 
the discrete sum over £ by an integral over £ or over the impact parameter 

Yl^ j d ^ k omj db, (23) 

the forward scattering amplitude can be written in terms of the impact parameter 



R 



f(0 = Q) = ^ I ( e *2(fc in -fc out )Kcos^ _ y 6d6; 



or in terms of the angle ip 



f(0 = 0) 



koutR 2 
2i 



^ _ e i2(k in -k out )Rcos^ dcos 2^,_ 



Then, the expression for the total cross section is obtained from Eq. 

w = cosip, we get 



Otot = Im 



Im 



Jo 



2nR 2 -i4iiR 2 / we l{ki »- k ™ t)2Rw dw 



Then, from integration by parts, we obtain 

e 



Otot = Im 
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i 2nR z - i 4nR z 



i 2ttR 2 - i AttR 2 
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i(k\ n fc ou t) 2R 

i (/cj n /c ut )2fi 
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1 e i(k in -k out )2Rw 

JO ^(^in — k out ) 2-R 
e i(fci n -fc ou t)2J? _ -y 



(24) 



(25) 

and setting 



(26) 
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i(k in -k ont )2R (2Rf(k in -k 



out J 



(27) 
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Fig. 9. (Color online) The extinction efficiency factor versus the phase shift Ar, for non- 
absorbing spherical particles with relative refractive indices n = 1.05, 1.33, and 2.00. The 
Ramsauer solution from Eq. (130]) for n = 1.05 is also given. The vertical scale applies only 
to the lowest curve, the others being successively shifted upward by 2. 

If we introduce the parameter Ar = 2R {k- in — k ont ), the analytic expression for the cross 
section becomes 



o- tot = 2nR 2 



sin Ar 
2 . - + 



sinA R /2 



(28) 



A R V A R /2 

In this expression, the total cross section approaches 2ttR 2 at high energies. This is the 
so called "black disk" approximation. Nevertheless, from a purely geometrical viewpoint, 
between two spheres, a collision occurs if the distance between their centres is less than ri+r 2 , 
where r\ and r 2 are their radii. So the geometrical cross section is equal to 7r(ri + r 2 ) 2 . In 
our case, the photon may be considered to be a particle with a diameter equal to its reduced 
wavelength A = \/2ir [20j. Consequently, Eq. (|28|) becomes 



2n[R + ^- 

47T 



2 sin A R ( suiAr/2 



(29) 



A R ' V Ar/2 

Finally, an expression for Ar as a function of the index of refraction n is needed. Ac- 
cording to the conservation of phase, it is k u 
X = 27r//c out 



nk mf . With this relation and the fact that 



Q 



ext,R 



irR 2 



2 1 



n 



1 



A 



R 



1 



sin Ar 
2 . + 



A 



suiAr/2 
Ar/2 



(30) 
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Fig. 10. (Color online) Total cross section as a function of the phase shift Ar = 2x(n — 1) for 
a refractive index fixed at (a) n = 1.05 and (b) n = 1.33. The simplified Ramsauer solution 

n 

from Eq. (I3TT) is in blue continuous line, the Mie prediction in red dashed line (using [161 ]) 
and the formula from Eq. (1281) in black thick line. The incident wavelength is fixed at 0.6 /jm. 

where Ar = 2R (k in — k ont ) = 2R {n — 1) A; out = 2x (n — 1). This result is derived using the 
approximation of scalar light, unlike the Mie solution. Therefore our formula applies only to 
scattering of unpolarised light. Note that the first factor of Eq. fl30|) explains the fact that 
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the amplitude of the large oscillations grows with n. The extinction efficiency factors, when 
plotted against Ar, show a universal shape (see Fig. [9]). Our model and the Mie predictions 
are in good agreement for the three refractive indices used. 
B. Ramsauer total cross section via the differential cross section 

The simplified Ramsauer formula is very convenient to use, in computing the phase function, 
as we have just seen, and also in calculating the total cross section. In this section, we 
compare the Mie prediction with the one obtained by integration of Eq. (ITS]) . From the 
differential cross section, the total cross section is given by 

2n it 

a tot = J d$ J dflsinfl^. (31) 
o o 

Fig. QUI compares, for two values of the index of refraction, the integrated Ramsauer phase 
function and the Mie theory result. The analytical solution from Eq. ff28]) is also plotted. 
Each curve is characterised by a succession of maxima and of minima. While there is 
generally a good agreement between the three curves, only the integration of the simplified 
Ramsauer solution PgR reproduces the Mie curve at small phase shift Ar, the so-called 
Rayleigh regime. The agreement improves when the refractive index decreases. In fact, at 
low phase shift, i.e. at small size parameter, the behaviour of Eq. ( 1251) is not proportional 
to A~ 4 as expected by the Rayleigh theory, but only to A~ 2 . Of the two methods, only the 
integration of the simplified Ramsauer solution is able to reproduce correctly the physics of 
Rayleigh regime. 

4. Application to the Henyey-Greenstein scattering phase function 

Monte Carlo simulations are usually used to describe light scattering in a medium. It is a 



very glo 
biology 



21] as well as in meteorology 



221 or 



problematic encountered in astrophysics 
23| . One of the critical points is to simulate the multiple scattering of photons in 
the medium. Also, the size distribution of particles or their shape are often unknown: the 
Mie phase function has to be evaluated for each particle. One of the most popular scattering 
phase functions is the Henyey-Greenstein (HG) function [24]. This function was first intro- 
duced by Henyey and Greenstein in 1941 to describe scattering processes in galaxies. This 
is a parameterisation usually used to reproduce scattering on objects large with respect to 
the incident wavelength. In the case where the backscattering can not be neglected, the HG 
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Fig. 11. (Color online) Henyey-Greenstein functions representing the scattering phase func- 
tion for different asymmetry parameters ^hg and backward factors /. The Rayleigh phase 
function, proportional to (1 + cos 2 #), is also plotted. 



function becomes a "Double HG" and is given by 



iW%HG,/) 



9hg 



47T 



;i + g nG 2 - 2shg cos Of 2 "2(1+ shg 2 
where 9 is the scattering angle, / represents the strength of the second component relative 
to the backward scattering peak (/ = meaning no additional component) and #hg is the 
so-called asymmetry parameter. It is defined as 



+ f- 



3 cos 2 e 



2^3/2 



(32) 



9kg = (cosf9) = / cos0P H G(0|0HG,/)27rsin0d0. 



(33) 



Henyey-Greenstein functions for different values of guc and / are plotted in Fig. [TTJ The 
asymmetry factor affects the strength of the forward peak and its width. Contrary to the 
Rayleigh phase function, Henyey-Greenstein phase functions depict a strong forward peak 
directly linked to the asymmetry parameter #hg- Using the Ramsauer approach developed 
in Sect. |2] and the resulting relation A6>fwhm = 2//c out i?, the mean radius of the particle 
size distribution can be estimated from the width of the forward peak, for a fixed incident 
wavelength. Fig. [T2] shows the relation between gnc values and mean particle radius for 
two different incoming wavelengths A. To take into account the overestimation of the mean 
radius explained in Sect. 12 B\ the values plotted have been reduced of 30%. In some cases, 



20 




0.4 0.6 OA 

Asymmetry parameter g H = < cos 6 > 

Fig. 12. (Color online) Relation between the asymmetry parameter of the Henyey-Greenstein 
function and the mean radius of the particle size distribution. The equivalence is plotted 
for two different incident wavelengths: 0.4 /iin (blue circles) and 0.8 /im (red squares). 

only the asymmetry parameter is measured by an experimental setup: these curves can be 
useful to estimate directly a mean particle size for the medium studied. If the first goal is to 
estimate the particle size distribution, a more detailed analysis seems to be needed. However, 
this method gives a first estimation of the particle size and an intuitive understanding of 
the asymmetry parameter $hg for the Henyey-Greenstein parameterisation. 

It is typical of what is needed in large ground-based astrophysics projects as cosmic ray 



shower detectors 25( or imaging atmospheric Cherenkov telescopes 26j. In these experi- 
ments, the atmosphere is used as a giant calorimeter and the associated systematic uncer- 
tainties have to be as small as possible 271 . Usually, atmospheric effects are divided into two 
categories: the molecular component [28], monitored using ground-based weather stations 
and meteorological radio soundings, and the aerosol component [29]. Atmospheric aerosols 
are small particles as dust or droplets in suspension having a typical size varying from a few 
nanometres to a few micrometres. Most of the atmospheric aerosols are present only in the 
first few kilometres above the ground. Unlike the molecular component, the aerosol popula- 
tion is highly variable in time and location, depending on the wind and weather conditions. 
Two main physical quantities have to be estimated to correct the effect of the aerosols on 
the number of photons detected by the telescopes: the aerosol attenuation length and the 
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aerosol scattering phase function. Some of these astrophysics projects have already pub 
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3jj, giving a first estimation of the 



lished an estimation of the asymmetry parameter 
mean aerosol size present in the low part of the atmosphere. 

5. Conclusion 

We have shown that within the experimental errors the Ramsauer method predicts an ex- 
tinction efficiency factor and a scattering phase function compatible with those given by the 
Mie solution. The large oscillatory behaviour for the extinction efficiency factor is under- 
stood as the consequence of the interference between the fraction of light going through and 
the fraction passing close to the sphere. 

Two Ramsauer solutions have been given for the scattering phase function. The form 
of the first one, Pr(9), without approximations, itself does not lead to obvious physical 
interpretation. The second, the simplified Ramsauer solution Psr(#), is an analytical formula 
from which characteristic behaviour can be predicted. However it cannot reproduce the 
Rayleigh regime for small size parameters. 

The Henyey-Greenstein function is a parameterisation intensively used in Monte Carlo 
simulations describing light scattering. Its key parameter is the asymmetry parameter ^hg- 
Using the Ramsauer approach, we proposed a basic relation between #hg an d the mean 
radius of particles scattering light. This approach is already employed and seems to be 
useful in large astrophysics projects. 
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